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SECTION  I 


INTRODUCTION 

Digital  controls  are  the  subject  of  expanding  interest  in  both  the 
academic  and  industrial  communities.  This  is  especially  true  since  the  advent 
of  the  economical  microprocessor  and  its  related  support  chips  which  have 
nade  digital  control  an  attractive  alternative  when  considering  a  control 
strategy.  For  the  past  few  years  new  control  systems  have  been  designed 
using  digital  instead  of  analog  controllers.  There  is  also  a  need  for 
converting  systems  which  were  designed  to  be  controlled  by  continuous  controllers 
into  systems  to  be  controlled  by  digital  controllers.  Conversion  of  existing 
control  systems  is  a  problem  of  considerable  interest  to  various  branches  of 
industry . 

A  popular  method  for  converting  existing  control  systems  into  digital 
control  systems  is  by  Tabak.*  This  method -converts  the  transfer  function  of 
the  continuous  controller  into  a  pulse-transfer  function  of  the  digital 
controller  using  prewarped  bilinear  transformation.  The  method  is  straight¬ 
forward  and  easy  to  use,  but  the  digital  control  system  obtained  by  this  method 
approximates  the  continuous  control  system  only  when  the  sampling  frequency  is 
sufficiently  high  as  compared  with  the  highest  frequency  of  the  continuous 
system.  This  is  because  the  tranfer  function  of  the  plant  and  the  feedback 
element  are  not  taken  into  consideration.  Thus  the  capabilities  of  the  digital 

controller  are  not  fully  utilized. 

2  3  4 

Kuo  et  al.,  Yackel  et  al.  and  Singh  et  al.  in  their  papers  proposed 
methods  for  converting  a  continuous  control  system  into  a  sampled-data  control 
system  by  matching  the  unit-step  response  of  the  systems  at  all  sampling  instants 
or  at  multiples  of  sampling  instants.  The  matching  requirements  are  satisfied 
by  changing  the  input  and  feedback  gains  of  the  system  instead  of  using  a 
digital  controller.  However,  such  matching  holds  only  for  unit  step  input.  Thus 
the  performance  of  the  sampled-data  system  obtained  by  their  method  approximates 
that  of  the  continuous  system  only  when  the  input  frequency  is  sufficiently  low 
in  comparison  with  the  sampling  frequency.  For  higher  input  frequencies,  the 
approximation  is  poor. 


Recently  the  author  of  this  report  and  Yeh^  proposed  a  computer-aided 
method  for  converting  an  existing  continuous  control  system  into  a  digital 
control  system  by  matching  the  discrete  frequency  response  of  the  discretly 
excited  system  to  the  frequency  response  of  the  continuous  model  with  a  mininura 
weighted  mean  square  error.  The  discrete  frequency  response  of  the  digital 
control  system  was  obtained  by  inserting  an  artificial  sampler  of  duration  T 

i  .  .ffi 

at  the  output  of  the  system  and  then  substituting  z  ■=  eJ  in  the  overall 
pulse- transfer  function  of  the  system.  Even  though  the  results  obtained  by 
this  method  were  superior  to  the  ones  obtained  by  Tabak's  and  other  methods, 
this  method  only  matched  the  response  of  the  systems  at  the  sampling  instants. 
The  purpose  of  this  study  is  to  improve  on  this  method  by  matching  the' 
continuous  frequency  response  of  the  discretly  excited  system  to  that  of  the 
continuous  model.  The  synthesis  of  the  digital  controller  is  carried  out  by 
minimizing  an  error  between  the  continuous  frequency  response  of  the  digital 
and  existing  systems  in  the  w-domain.  The  method  uses  the  complex  curve 
fitting  technique  of  Levy.**  As  a  result  the  formulas  for  computing  the 
parameters  of  the  digital  controller  are  obtained.  This  method  is  then  extended 
to  the  digital  redesign  of  an  existing  multiloop  continuous  control  system. 

The  design  algorithms  developed  here  are  based  on  the  transfer  function  of  the 
plant  and  the  closed-loop  transfer  function  of  different  loops  of  the  continuous 
system.  \s  a  particular  application,  the  digitalization  of  a  flight  control 
system  is  discussed.  The  system  under  consideration  is  the  flight  controller 
for  the  longitudinal  YF-16  aircraft.  The  results  obtained  by  this  method  are 
compared  with  those  obtained  by  the  previous  methods. 

SECTION  II 

STATEMENT  OF  THE  PROBLEM 

The  block  diagram  of  a  general  multiloop  continuous-data  control  system  may 
be  drawn  as  shown  in  Figure  1 


CONTINUOUS 

CONTROLLER 


Gc3(  ®  ) 


Figure  1,  Existing  multiloop  continuous-data  control  system 


The  continuous  model  shown  in  Figure  1  can  be  digitalized  by  inserting 
a  sampler  at  the  error  input  e(t)  and  replacing  the  continuous  compensators 
Gc^(s)  by  digital  controller  D^(z)  as  shown  in  Figure  2. 


Figure  2.  Multiloop  digital  control  system 
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where  a’s  and  b’s  are  unknown  constants.  It  is  required  to  design  digital 
controllers,  i.e.,  to  find  the  unknown  constants  of  different  controllers 
so  that  the  steady  state  output  of  the  digital  control  system  follows  the 
desired  output  of  the  continuous  model  for  all  sinusoidal  inputs  within  the 
frequency  range.  This  is  done  by  matching  the  continuous  frequency  responses 
of  different  loops  of  the  digital  control  system  to  the  corresponding  ideal 
frequency  responses  starting  from  the  innermost  loop. 


SECTION  III 


DESIGN  OF  DIGITAL  CONTROLLERS  VIA  CONTINUOUS  FREQUENCY  MATCHING 

1.  Digital  Controller  in  the  Feedback  Loop 

The  block  diagram  of  the  inner  loop  of  the  digital  control  system  shown 
in  Fig.  2  is  redrawn  as  shown  in  Figure  3. 


PLANT 


Figure  3.  Inner  loop  of  digital  control  system  shown  in  Figure  2 


The  overall  transfer  function  of  the  continuous  and  digital  systems  (inner 
loop)  can  be  written  from  Figures  1  and  3  respectively  as 


h 


(2) 


Cml<8>  _  ,  ,  C1(S) 

R^s)  “  'lW  “  l-KJ1(s)Grl(s) 


Cdl(8> 


Gdl<s) 


GhOGl(s> 

i+dI(s)Vit<s) 


(3) 


where  G  _(s),  G^q(s)  and  G^(s)  are  the  transfer  functions  of  the  continuous 

controller,  zero-order  hold  and  plant  respectively.  The  superscript  notation 

T 

denotes  the  sampling  rate.  For  example  indicates  that  the  signal  R^  is 
sampled  at  samples/second. 

It  is  required  to  design  a  digital  controller,  D^(z),  by  matching  the 

continuous  frequency  response  of  the  digital  control  system  shown  in  Figure 

3  to  the  corresponding  ideal  frequency  response  of  the  continuous  model 

shown  in  Figure  1.  The  frequency  response  of  the  continuous  model  can  be 

obtained  by  substituting  s  =  ju  in  equation  (2)  and  varying  u>  from  0  to  ». 

But  because  of  the  hybrid  nature  of  equation  (3),  a  common  way  of  finding  the 

frequency  response  of  the  digital  control  system  is  to  put  an  artificial  sampler 

of  duration  T  seconds  at  the  output  of  the  system,  substituting  z  =  in 

the  overall  z-transfer  function  and  then  varying  u  from  0  to  ir/T.  We  will 

call  this  discrete  frequency  response.  Recently  the  Author^  proposed  a  method 

for  converting  an  existing  multiloop  continuous  control  system  into  a  digital 

control  system  by  matching  the  discrete  frequency  responses  of  the  digital 

control  system  to  the  frequency  responses  of  the  continuous  model.  In  this 

method  we  only  matched  the  response  at  the  sampling  instants.  It  is  this 

observation  which  prompted  the  development  of  the  method  in  this  report 

which  matches  the  continuous  frequency  responses  of  the  discretly  excited 

system  shown  in  Figure  1  to  those  of  the  continuous  system  shown  in  Figure  2. 

s  t 

Rewriting  equation  (3)  with  z  =  e  ,  we  get 


Gdl(s) 


GhOGl(s) 


(A) 


1+Di(z)GhoGi(z) 
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Then,  the  continuous  frequency  response  of  this  discretly  excited  system  is 

8 

given  by  (Didaleusky  and  Whitback  ) 


where 


A  +1B 
n  J  n 


f  Viw 


s-Jo,  J+D(z)ChOGl(z) 
n 


A  sT 
z  *»  e 


z  =  bT 


2nn 


n  **  0,1, 2,3, 


(6) 


(7) 


For  n=0,  we  get  the  fundamental  component  of  the  output  response  (i.e.  "'^(s) 
for  a  sampled  sinusoidal  input  R(z)).  For  n^O,  we  get  the  output  res  se  for 
the  individual  alias  frequency  components.  For  example,  for  n=l 


u 


1 


b  +  a) 

s 


(8) 


and  we  get  the  frequency  response  in  the  output  for  the  1st  alias  component. 

At  this  point  it  is  important  to  recognize  that  the  actual  continuous 
response,  C^(s)  ,  is  not  sinusoidal.  This  becomes  clear  in  the  time 
domain  response  for  C^(s)  for  a  sampled  input  R^(z)  (Whitback  and 
Hoffman^).  It  is  apparant  that  the  best  we  can  do  is  match  the  fundamental 
component  of  the  output  response  C^(s)  .  An  alternate  w-domain  expression 
for  equation  (6)  is  given  by 


A  +1B 
n  J  n 


f  GhOGl(s) 


S=jU) 


ga<"> 


n 


,  .2 

'w=3Y 


tan 


w  T 
n 


(9) 


where 


Ga(w) 


_ 1 

l+D(z)Gj^G1(z) 


(10) 
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Equations  (6)  and  (9)  provide  the  same  continuous  frequency  response  of  the 
dlscretly  excited  system  but  the  advantage  of  using  equation  (9)  is  that  it 
improves  numerical  round  off  problem  in  the  implementation  algorithms.  Equation 
(9)  can  be  written  as 


Cdi(w)  **  — W  - - - - 

1  +  ^  l+D(z)Gh0G1(z) 


GjCs) 


wT 

1  +  “ 

1  -  wT 
2 


.2  -1  wT 

|s=Jy  tan  2J 


.  Let  the  transfer  function  of  the  digital  controller,  D^,  be  written  in  the  w-domain 


D1(w) 


TO,  TO— 1. 

a,  w  4a-  .w  4  ' 

1  ,m _ 1  ,m-l _ 

.  n.,  n-1, 

b,  w  4b.  .w  4- ■ 

l,n  l,n-l 


where  a1  's  and  bn  's  are  constants  to  be  determined.  The  continuous  frequency 

1,1  1,1 

responses  of  the  digital  and  continuous  system  can  now  be  obtained  by  substituting 
2  -l  yT 

w  «  jy  an<*  s  =  j—  tan  A—  in  equations  (11)  and  (2)  respectively 


VJr)  -  Vs)  ,  _1T  -  a1+JBi 

r=M —  tan 


\s~i~  tan 


where 


Gdi(3Y)  - 


<Xll^Yll) 

Pi+jyQi 


1  +  14L14jyM1  (X21+jY21) 


P1  "  al,0-ai,2Y2+al,4Y^~al  4  •••• 

2  A  6 
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1  1,2  1,4*  1,6 


"  bl  i“bi  ^  +bi  cY  -b*  -7Y6+ 
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xn+jYn  ' 


jb_  ci(s) 

14.<U-  * 


X21+JY21  *  GhOGl<2) 
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T 
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T 
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For  notation  convenience,  let  the  numerator  and  denominator  of  G ,, (1y)  be 

ui 

denoted  by  N(y)  and  M(y)  respectively,  i.e. 

s' 

Gdi(3Y>-Rw 

In  terms  of  real  and  imaginary  parts,  N(y)  and  M(y)  can  be  written  as 

N(Y>  -  4«+j0 
M(y)  *=  Cf+jT 


where 


*  "  (l+L1)XirYM1Yu  <2< 

0  -  (1+Lpy^+YM^j^  C2I 

«  '  <2< 

The  error  between  G^^Cjy)  and  F^(.ly)  is  defined  as 

e(Y)  “  F1(jy)-Gdl(jY)  (2£ 

An  optimal  set  of  a’s  and  b’s  can  therefore  be  optained  by  minimizing  the  mean 
square  of  the  magnitude  of  e(y)  Let  the  error  function  to  be  minimized  be 
defined  as 


« 


(29) 


J> 


Uy>  - 


where  (Y2-Y^)  *s  tbe  bandwidth  in  the  w-dora  'n  on  which  the  matching  is  required 
This  error  function  is  to  be  minimized  by  adjusting  the  unknown  coefficients 
*1,1  and  bi  i*  However»  due  to  the  presence  of  M(y)  in  the  denominator,  the 
equations  obtained  by  differentiating  E  with  respect  to  a.  and  b  are  nonlinear. 

This  nonlinearity  may  be  removed  by  modifying  the  error  function  (29)  as 


r2 

iFi 

«  V 


(Jy)M(y)  -  N(y) |  dy 


This  integral  differs  from  that  in  equation  (29)  by  including  |M(y)|2  in  the 
Integrand  as  a  weighting  factor.  This  modification  is  not  new.  It  was 
suggested  by  Kalman^ as  an  identification  technique,  was  used  by  Levy  b  in 
complex  curve  fitting,  Rao  and  Lamba^  for  simplifying  linear  systems  and  by 
the  author  and  Yeti'*  for  discretizing  continuous  systems. 

The  integrand  in  equation  (3)  is  a  complex  function  and  may  be  seperated  into 
real  and  imaginary  part  as 

|F1<jY)M(Y)  -  N(y)j  *  A(y)  +  jm(y)  (31) 

where 

Uy)  *=  AJLo-B1T-(f  (32) 

m(y)  *=  AjT+BjO-0  (33) 

By  substituting  equations  (31) —(33)  in  (30),  we  get' 


£(A1a-B1T-<|>)  ^(AjT+BjO-G)  2]  dy 


Differentiating  Effll  with  respect  to  a ^  1's  and  bx  and  equating  the  resulting 

equations  to  zero  represents  (mfn+1)  linear  equations  in  (m+n+1)  unknowns 

„  to  a,  and  b,  to  b,  .  These  equations  can  be  written  in  the  matrix 
1,0  l,m  1,1  l,n 

form  as  given  in  equation  (35) 
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0  -T4  0  T6  0  -Tj- 


Vi  W2  "V3  'W4  V5  V 

“U2  v3  ~V5  -u6  V 

-V3  -°4  V5  °6  *V7  ‘V 

U4  *V5  -°6  V7  °8  ~V9* 


V1  *°2  “V3  °4  V5  ”D6 . Uo 

u2  v3  -V  -v5  u6  -v, .  vi 

U4  V5  -06  -V7  »*•••••  u2 

"V5  U6  V7  -°8  _V8 . 


w2  o  -v4  0  Vfc  0 

o  w4  o  -w6  o  w8 

-w4  o  w6  0  -v8  0 

o  w6  0  W8  o  -Wl0 


where 


(AK)  (4l+Y2l)^ 


1+BJ)X21-A1  (X11X21+Y11Y2^  +B1  (Vai'Vll)]  Y”dh 


°h  ■  h 


[(AX)Y2rAl(: 

*  v 


l)-Bl(X2 


XllY2l“X21Yll)_BllXllX21+YllY2l)  Y  dh 


wh  •  2  [fi +BJ)+  (Xn+Yii)  -2  (aixii+biyu)]y>,'Jy 

s 
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Equations  (35)— ( 39)  provide  a  design  algorithm  which  has  been  programmed  on 
a  digital  computer  (IBM  370/3031)  and  is  available  in  the  form  of  a 
subroutine  written  in  FORTRAN  IV.  Note  that  the  order  of  the  digital 
controller  (integer  m  and  n)  is  arbitrary.  Therefore,  using  this  algorithm, 
one  may  try  out  several  digital  controllers  of  different  orders  and  select 
the  best  design  by  considering  the  tradeoff  between  cost  (complexity)  and 
performance.  The  pulse-transfer  function  of  the  digital  controller  can 
now  be  obtained  as 

D1(z)«D1(w)|  2  z_x  (40) 

w  =  T  z+1 

2.  Digital  Controller,  D2(z),  in  the  Feedforward  path  of  Outer  Loop 

The  block  diagram  shown  in  Figure  4  is  obtained  from  Figure  1  after  inserting 
a  sampler  at  the  error  input  e(t)  ,  replacing  continuous  compenstor  ^^(s) 
by  the  already  designed  controller  D^(z)  and  replacing  the  controller  G  ^(s) 
by  the  yet  to  be  designed  digital  controller  D2(z). 


R(l) 

R(s) 


DIGITAL 

CONTROLLER 

D2(z)  — ' / 


PLANT 


DIGITAL 

CONTROLLER 

—  D,(z) 


CONTINUOUS 

CONTROLLER 

—  GC3(s)  -**- 


C(t) 

C(s) 


Figure  4.  Digital  control  system  with  unknown  digital  controller 
in  the  feedforward  outer  loop. 
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The  overall  transfer  functions  (outer  loop)  of  the  digital  control 
system  as  shown  in  Figure  4  and  the  continuous  control  system  shown  in 
Figure  1  can  be  written  in  w-domain  as 


F2(w) 


Gc2^F1^G2^ 

1+Gc3(s)Gc2(s)F1(s)G2(s)  ,  -1  wT 

8“jj  t®n  2j 


Gd2(w) 


,  V2>  ltChOCl(z)Bl(2> 

^  1+Vz)  Sflpig .  „ 

1+GjiqG^(z)D^(z)  1+- j 

2 


(s> 


.2  -lwT 
|s-j-  tan  2j 


where 


GhOGlG2<2>  *  z[GhOGlG2(s)] 

'  2  [GhOGlG2Gc3(s)] 


Let  the  transfer  function  of  the  digital  controller,  D2(w),  be  written 


in  the  w-domain  as 


m,  m-1. 

a„  w  +a„  .w  + 

_  „  x  2.m  2, m-1 

D  (w)  »  — 1 — - rrr- 

2  b„  wn+b-  .Wn  K 

2,n  2,n-l 


+  Vl^.O 

+  b2(1vfl 


,  m<n 


where  a2  ^'s  and  b2  ^'s  are  constants  to  be  determined.  The  continuous  frequency 

responses  ot  the  digital  and  continuous  system  are  obtained  from  equations 
(41)  and  (42)  by  substituting  w  =  jy  and  are  given  by 

G.0(jto)F1(ja))C0(ju))  I 


l+G-UuOC^OfaOF,  (ju>)G„(ju>) 


where 


l+L24jYM2  y  12  J  12J 

Gd2(jY)  “  P2+jYQ2  /  \ 

1-ll+L2+jYM2  p2+1Y22j 


P2  “  a2,0"a2,2Y2+a2,4Y4"a2,6Y6+  **“ 

2  4  6  • 

Q2  a2,l“a2,3Y  +a2,5Y  ~a2,7Y  +  **** 


L2  "  °"b2,2Y2+b2,4Y4"b2,6Y6+ 


H2  "  b2,rb2,3Y2+b2,5Y4-b2,7Y6+  *'** 


X12+JY12 


l+jYj  1+Gh0Gl(z)Dl(z) 


(jtu) 

wT  jw 
2 


2  -1  yT 

tan  -J- 


X22+JY22  * 


GhOGlG2Gc3(z) 

l+Gh0GjL(z)D1(z) 


T 

i+jrj 

■  t" 
1-JT2 


The  error  function  between  G^Cjy)  and  ( j Y)  which  is  to  be  minimized  is 

given  by 


h°- 


B2t-*)  +  (A2t+B2o-0)  JdY 


j-0)2]< 


where 


*  "  P2X12+yQ2Y12 
0  "  P2Y12^2X12 


o  -  1+L2+P2X22 


T  -  ym2+p2y22+qx22 


Differentiating  E  with  respect  to  a_  ,'s  and  b_  ,'s  and  equating  the 

D*  f  1  4)1 

resulting  equations  to  zero  represents  (o+n+1)  linear  equations  in  (nrfn+1) 

unknowns  a„  n  to  a.  and  b,  1  to  b  .  Then  equations  can  be  written  in  the 
2,0  2, in 

matrix  form  as  given  in  equation  (59) 


0  -T, 


-T2  0 
0  -T, 


0  "T, 


T*  0 


V1  W2  *V3  “U4  V5  V 

-u2  v3  uA  -v5  -u6  v7- 


°4  ~V5  -°6  V7  U8  ~V9* 


*1  -*2 


-°4  -V5 


-v3  v5  *»6 

-°4  -V5  U6  -V7 
V5  "U6  -V7  U8 
C6  V7  -°8  ‘V8 


u2  0  -w4  0  w6  0 

o  «4  o  -w6  o  w8 


-v3  -««  V5  U6  -V7  -Ufi . -VA  0  V6  0  -W#  0 


o  vfc  o  w8  o  -w10 


where 


ri(AH(x»+ 

yi 

\  -  P(AH)X22-(A 


X22+Yy  +  (X12+Y12, 


>)+2B2  (! 


,X12Y22-X22Y12j 


.)-  2A2(: 


lX12X22Y12Y22 )  Y  dY 


2X12+BY12  Y  dY 


Vh  -  f(A2+B2)Y22+(B2X12-A212)- 

Vh  ■  f2  (A2+B22)y1,<!Y 


Equations  (59)- (63)  provide  a  design  algorithm  which  has  been  programmed  on  a 
digital  computer,  IBM  (370/3031)  and  is  available  in  the  form  of  a  sub¬ 
routine  written  in  FORTRAN  IV.  The  pulse  transfer  function  of  the  digital 
controller,  D^(z),  can  now  be  obtained  as 

D2(z)-D2(w)|  2  Z_1  (64) 

wf  z+1 

3.  Digital  Controller,  D^(z) ,  in  the  Feedback  Path  of  Outer  Loop 

The  block  diagram  of  the  digital  control  system  after  replacing  the 
continuous  compensators  Gc^(s)  and  G^te)  with  already  designed  digital 
controllers  and  after  replacing  Gc^(s)  with  yet  to  be  designed  digital, 
controller  D^(z)  is  shown  in  Figure  5. 


PLANT 


C(l) 


CU) 


Figure  5.  Digital  Control  System  with  unknown  digital  controller  in  the 
feedback  loop. 
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The  overall  transfer  function  t  the  digital  control  system  shown  In 
Figure  5  can  be  written  in  the  w-donwin  as 


Ca3(w) 


n2(z) 

wQGiWV.) 

?  D  (z)G  G  G  (z) 

2  1+D,(z)  -2--- K°-L-2 - 

1+GhOGl(z)Dl(2) 


G1G2(s) 


's«j|  tan-1 


Let  the  transfer  function  of  the  digital  controller  D^(w)  be  written  in 
the  w-domain  as 


Vw> 


in.  m—1 

a3.mw  +a3.m-lw  + 

N  nWn+b-3  n 


a3. 1^3.0 

+  V^.o 


,  m<n 


where  a^  ^'s  and  b^  ^'s  are  constants  to  be  determined.  Since  the  digital  controller 

D^(z)  is  located  in  the  feedback  path,  these  constants  can  be  determined  using 
the  design  algorithm  developed  for  D^(z),  equations  (36)-(41),  where 


B1  **  B2 


xii+jYn  -  T~wT 


Dz(2) 


G1G2(s)| 


1+=J  l+Gh0G1(z)D1(z) 


,wT  s 


I  • 2  -  -I  YT 
1  s~j-  tan 


X21+JY21  •“ 


D2(2>GhOGlG2(z) 

l^Gl(z)Dl(z) 


,  ,wT 


Finally,  we  would  like  to  point  out  here  that  using  the  design  algorithm 
developed  for  the  digital  controllers  in  the  feedforward  and  feedback  path, 
we  can  convert  any  existing  ir.ultiloop  continuous  control  system  into  a 
multiloop  digital  control  system. 
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SECTION  IV 


NUMERICAL  EXAMPLES 


1.  Numerical  Example  1 

The  first  example  considered  in  this  report  is  to  convert  an  existing 


single  loop  continuous  control  system  shown  in  Figure  6  into  a  digital 
control  system  shown  in  Figure  7. 


Figure  6.  Existing  single-loop  continuous-data  control  system  for 
Example  1. 


Figure  7.  Redesigned  digital  control  system  for  Example  1. 

The  transfer  function  of  the  continuous  controller,  plant  and  the  feedback 
element  is  given  by 


„  /  x  1+0. 416s 

°C(B)  '  1+0. 135s 

(71) 

x  10 

C(s)  “  s(s+l) 

(72) 

H(s)  «  1 

(73) 
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Using  equations  (71)-(73),  the  overall  transfer  function  of  the 
continuous  model  can  be  written  as 


F(s) 


Gc(s)G(s) 

l+Gc(s)G(s)H(s) 


29.93s  +  71.942 _ 

s3+8.194s2+37.12s+71.942 


(74) 


The  magnitude  of  the  frequency  response  of  the  continuous  model  is  obtained 
using  equation  (74)  and  is  shown  in  Figure  8. 


Figure  8.  Magnitude  of  the  frequency  response  of  the  continuous 
system  of  example  1. 

Based  on  the  frequency  response,  the  sampling  peribd  T  is  selected  to  be  0.15s 
From  equations  (72)-(73), 

10 (l-e"sT  ) 
s2(s+l) 


GhGH(s) 


(7 


The  z-transform  associated  with  this  transfer  function  is  given  by 


G  gH(z)  =  0>10708z+0. 101858 
h  z2-l . 8607  z+0 .8607 

1+jy? 

Substituting  z  =  - —  in  equation  (76)  we  get 


9 .981-. 73w-.0019w 


chGH(3r)  -  x22+jy2  = 


w=jy 


From  Figure  7  and  equation  (72)  ,  we  get 

J.Y.. 

T 

1+J*2 


*12^12  - 


lY _ IQ 

l+jrf  s2(s+i) 


.2  -1  YX 

8  "  tan 


(76) 


(77) 


C78) 


8  .  jf  tan'1  ^ 
JT  2 


Also 


.  .  29.93s+71.942 

A  +jB  =  — - 

sJ+8.194s+37.12s+71.942 


s=j|  tan-1 


(79) 


Equations  (77)— (79)  are  the  only  information  required  to  design  a  digital 
controller  for  the  digital  control  system  shown  in  Figure  7.  Feeding  these 
equations  to  the  computer  program  written  for  equations  (59)-(63),  the 
transfer  function  of  the  first-order  digital  controller  given  by  the 
algorithm  in  the  w-domain  is 


D(w) 


.59713  +  . 3636w 
l+.13053w 


(80) 


2  z— 1 

Substituting  w  =  —  in  equation  (80)  gives  the  pulse-transfer  function  of  the 
digital  controller 


D(z) 


1.9865  -  1.55074 
z-. 27108 


(81) 


The  frequency  responses  of  the  existing  continuous  system  and  the  digital 
control  system  with  the  digital  controller  given  by  equation  (81)  are  shown 
in  Figures  9-11.  It  can  be  seen  that  both  the  magnitude  and  phase  angle  of 


the  digital  system  matches  closely  with  that  of  the  continuous  model.  If 
the  designer  Is  not  satisfied  with  the  results  and  wants  to  improve  the 
performance,  the  only  change  which  need  to  be  made  is  to  increase  the 
order  of  the  digital  controller.  The  algorithm  goes  through  the  standard 
steps  given  above  and  finds  the  parameters  of  the  new  digital  controller. 
Thus  the  performance  may  be  improved  at  the  expense  of  increasing  complexity. 

2.  Comparison  with  the  existing  methods 

The  digital  controller  for  the  digital  control  system  discussed  in  the 

12 

numerical  example  is  redesigned  using  the  Tustin  transform  (bilinear 
transformation);  methods  presented  by  Tabak^  (prewarped  bilinear 
transformation)  and  by  the  author  of  the  paper  and  Yeh^  (discrete  frequency 
matching).  Using  these  methods,  the  continuous  controller  given  in 
equation  (71)  is  converted  into  digital  controllers  and  are  given  by 


Tustin  transform: 


Prewarped  bilinear  transform: 


2.294z-1.5936 

Z-.299 


w  ^  2.42965z-l .68046 

»<*>  ■  - 

Using  the  method  presented  by  the  author  and  Yeh: 


first-order: 


3.6787z-2.31256 
z+0. 383656 


second  order : 


D(z)  =  9,6.2  - 3 . 65z-f°  ■  56_  (85) 

z  -0.34z-0.384 

The  frequency  responses  of  the  digital  control  system  with  digital 
controllers  given  by  equations  (82)-(85)  are  shown  in  Figures  9-11.  Comparison 
of  the  results  in  these  figures  demonstrate  that  the  digital  controller 
obtained  by  the  present  method  gives  better  results.  The  unit  step  responses 
of  the  continuous  model  and  digital  control  systems  obtained  by  different 
methods  arc  shown  in  Figure  12.  It  can  be  seen  from  this  figure  that  the 
present  method  gives  better  results. 
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Nyquist  plot  of  the  frequency  responses  of  the  continuous  modei  and  the  digital 
control  system  designed  by  different  methods  (example  1). 


Magnitude 


i  I  i  I  I  r 


o 

o 

o 

o 

o 

o 

o 

o 

id 

o 

id 

o 

id 

o 

CD 

co 

CO 

CM 

N- 

• 

• 

T” 

1 

i 

CM 

• 

CM 

a 

seojBaa  ui  8|6uv  esei^d 


i 


23 


Figure  11.  Phase  plot  of  the  frequency  responses  of  the  continuous  model  and  the  digital 
control  system  designed  by  different  methods  (example  1) . 


—  Continuous  Model 

First-Order  Controller  (Present  Method) 

-  First  Order  Controller 

(Prewarped  Bilinear  Transform) 

-  First-Order  Controller 

(Tustin  Transform) 


_J _ i _ l _ 

1 .0  2.0  3.0 

Time  in  seconds 


Unit  step  responses  of  the  continuous  model  nnd  the  dipitnl 
control  system  designed  by  different,  methods  (example  1) . 


3.  Numerical  Example  2 

The  digitalization  techniques,  described4  in  section  III  was  applied 
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to  an  analog  flight  controller  for  the  longitudinal  YF-16  as  shown  in 
Figure  13.  A  listing  of  transfer  functions  of  the  aircraft  flight  at 
30,000  ft  at  M  *  0.6  can  be  found  in  Appendix  A. 

The  problem  is  to  digitalize,  i.e.,  to  replace  the  continuous-time 
controllers  by  digital  controllers  while  preserving  the  dynamic  properties 
of  the  whole  system.  The  block  diagram  of  the  digitally  controlled  YF-16 
obtained  by  inserting  a  sampler  at  the  input  of  the  actuator  and  replacing 
continuous  controllers  by  digital  controllers  is  shown  in  Figure  14. 


CONTINUOUS 

CONTROLLER 


G,(s) 

MODEL  OF^ 
ACTUATOR  I 

I  20  I  6.  I 


G*(s) 


G3(s) 


I - - 


CONTINUOUS  G  , 
CONTROLLER  cl 


CONTINUOUS  CONTROLLER 
GCj(  s ) 


CONTINUOUS  CONTROLLER 
Gc«(s) 

_  »3(s*5)  _  15  1 

s+15  s+15 


Figure  13,  Analog  YF-16  controller  (M=0.6,  h=30,000  ft) 


Dj(z)  Controller:  From  Figure  13,  the  transfer  function  of  the 
continuous  model  (innermost  loop)  for  the  design  of  digital  controller  D^(z) 


can  be  written  ns 
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c 


0,(2) 


Z.O.H.  II  20 
Gko(s)  s*20 


DIGITAL 

CONTROLLER 


_ X 

r\.t 

-  - 

J  I 

DIGITAL 

CONTROLLER 

) 

T 

D3(2) 


DIGITAL 

CONTROLLER 


Figure  14.  Digital  YF-16  controller 


(s)  l+Gcl(s)G1(s) 


1.45S4+162 . 1S3+1478 . S2+24 . 53S+8 . 3 


S  +30.95S-'+226.2S  +127.7S  +270.5S‘'+3.75S+1.31 
In  view  of  the  innermost  loop  of  Figure  14  and  equation  (19),  we  have 


X  +1Y  -  ^ - —  — (s) 

U  U  -4  S+2°  *•  .  -  4  If 


n  \  l~C~St  20  O  ,  , 

GhOGl^S  “  S  6+20  6  (  J 

e 

The  z-transform  associated  with  this  transfer  function  is  given  by 


eh0cx<*) 


.0022?*-.  00034s3-.  0057z2+.0037r.+  .  000175 
z5-4 . 41//*+7 . 68z  3-6 . 54  z 2+2 . 7  Is-  .4326 


2i 


1 


Substituting  Z 


T 

1+YjT 

- —  in  equation  (89),  we  get 

l-TJ 


Vl<lT)  *  *21+JY21 


.0029w5- . 05w*-l . 04w3+l AO . lw2+2 ■ 24w- . 79 
w5+19 .9w4+15 . 7w3-44 . 3w2- . 74w- .27 


'w-jy 


Feeding  equations  (86),  (87)  and  (90)  to  the  computer  program  written 
for  the  algorithm  implementing  equations  (36)- (40),  the  transfer  function 
of  the  first-order  digital  controller  in  the  w-domain  is  given  by 


Dx(w) 


0.501+.0109w 

l+.1007w 


2  z-1 

substituting  w  «  —  — -r-  in  equation  (91),  we  get 
T  Z+i. 


^  •  .1732z- .0072  . 

Dl(z)  *  Z-.668 

The  frequency  response  of  the  innermost  loop  of  the  digital  control 
system  shown  in  Figure  14  with  the  digital  controller  given  by  Equation  (92) 
is  shown  in  Figures  15,  16,  17.  To  compare  the  result,  digital  controllers 
are  also  obtained  by  the  other  methods  and  their  pulse-transfer  function 
are  given  by 


Tustin  transform: 


Dx(z) 


♦083(z+l) 

Z-.667 


Discrete  frequency  matching  method: 


D1(z) 


. 3659 (z-. 79 74) 
Z-.8553 


The  frequency  responses  obtained  using  the  digital  controllers  given  by 
equations  (93)  and  (94)  are  also  plotted  on  Figures  15,  16,  17.  Comparison  of 
these  results  demonstrate  that  the  digital  controller  obtained  by  the  present 
method  gives  better  results. 

D^(z)  Controller:  From  Figure  13,  the  transfer  function  of  the  continuous 
model  for  the  design  of  digital  controller  D2(z)  can  be  written  as 


F^(s) 


q(s)  Gc2(s)F1(s)G?(s) 

RjTs)  "  1+Gc  3  (  s  )  G  2  (  a  )  (  s )  C  2  (  s7 


_  lftOsS-.5r>8x104S5-4-.S79xl05sV.?Uxl06S3-K2.Vtx1(rS2+.7  3vl05S-+629.9 
S8447S/-l923S6+.86xlO/,S5+.35xloV,-»  .75xl(T53+.33xl0 V+329  1S+19.6 


(95) 


Continuous  Model 

Present  Method 

Discrete  Frequency 
Matching  Method 

Tustin  Transform 


-16.0  -12.0  -8.0  -4.0  0 

-5.0* 


4.0  8.0  12.0  Re 


Figure  15.  Nyquist  plot  of  the  frequency  responses  of  the  innermost  loop  of 

the  continuous  systcm  and  the  corresponding  digital  control  system 
designed  by  different  methods,  (example  2). 


Frequency  in  rad./sec. 


Figure.  16.  Magnitude  plot  of  the  frequency  responses  of  the  innermost  loop  of 
the  continuous  system  and  the  corresponding  digital  control  system 
designed  by  different  methods  (example  2). 


Phase  angle  In  degrees 


In  view  of  Figure  14  and  equation  (53),  we  can  write 


y  -My  b  —lX _ 

*12  3  12  ...  T 


1+jYj  l+Gh0GL(w)D1(w) 


T 

l+JYI 

T 

l-jJrf 


_20_  __q 

s+20  6e 


(s) 


(96) 


.2  _  -1  yT 

8"Jt  tan  2 


-ST 


G.  G. G_G  _(s)  -  -  -  24S_(S+5) 

no  12  c3'  J  S  (S+1HS+ 


S  .  (S+l) (S+15) (S+20)  6e 


(97) 


From  the  z-transform  of  this  transfer  function  and  equations  (90)  and 
(91),  we  get 


X22+jY22 


_ 1 _ 

1+Vi(w)Di(w) 


. 00003w6- . 004w5- . 36w VlO . 5w3+69  0 . 7w2+ll . lw+3 . 9 
w6+29 . 8w5+212 . lw4+110 . 6w3-436 . 7w2-7 . 6w-2 . 6 


k-jy 

(98) 


Feeding  equations  (95),  (96)  and  (98)  to  the  computer  program  implementing 
equations  (60)-(64),  the  transfer  function  of  the  digital  controller  D2  in 
the  w-domain  is  given  by. 


D2(w)vb  1.29  + 


(99) 


Substituting  w 


2  z-1 
T  z+1 


in  equation  (99) ,  we  get 


»2(z) 


1.399z-1.181 

(z-1) 


(100) 


The  pulse-transfer  functions  of  D2(z)  obtained  by  the  existing  methods  are 
given  by 


Tustin  transform: 


_  ,  ,  1.342(z-.8182) 

d2<z) - (i=i) 


Discrete  frequency  matching  method: 


D2(z) 


1.3837(z~. 80235) 
(z-1) 


(101) 


(102) 


The  frequency  responses  of  the  middle  loop  of  the  continuous  system  shown 
in  Figure  13  and  the  digital  control  system  shown  in  Figure  14  with  digital 
controllers  D2(z)  given  by  equations  (100)-(102)  are  shown  in  Figures  18,  19 
and  70.  Comparing  these  results  reveals  that  the  digital  controller  obtained 
by  the  present  method  gives  better  results. 
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Figure  20. 


Phase  plot  of  the  frequency  responses  of  the  middle  loop  of  the 
continuous  system  and  the  corresponding  digital  control  system 
designed  by  different  methods  (example  2). 


Dj(z)Controller 


From  Figure  14  and  equations  (70)  and  (71),  we  can  write 


xu+1Yu 


lx  Vu) 


1+->^2  1+ChOGl(u,Vu> 


-3— (s) 

s+20  6e'-s; 


w=jY 


.2  -1  T 

s=j-tan  yj 


(103) 


and 


X21+JY21  " 


I>2(w) 


- . 14w*+7 . 3w3+3 . 9  7w2+ . 034w 


1+GhOGl (w) °1 W  + • 95w ~2 • 33w  “ • 038w~ • 014 


w=jy 


(104) 


Feeding  equations  (95),  and  (103)  and  (104)  to  the  computer  program 
implementing  equations  (36)-(40),  the  transfer  function  of  the  digital 
controller  D_(w)  is  given  by 

J  » 

D,(w)  -  ^A1+096w-)— 

3  l+1.07w+.106w2 


(105) 


2  z-1 


substituting  w  =  —  in  equation  (105),  we  get 


D3(z) 


.8151z  -1.502z+.687 
z2-1.652z+.6646 


(106) 


The  pulse  transfer  functions  of  D^(z)  obtained  by  the  previous  methods 
are  given  by 

Tustin  transform: 


^  ,  0.995z  -1.81z+.814 

D  (z)  =  - - - - 

z  -1.5z+.517 


(107) 


Discrete  frequency  matching  method: 

„  ,  x  1.5433(z-. 8512) (z-1) 
D3(z)  "  (z~.3876)  (z-~9tT:>2) 


(108) 


The  frequency  responses  of  the  middle  loop  of  the  continuous  system  shown 
in  Figure  13  and  the  digital  control  system  shown  in  Figure  14  with  digital 
controller  D^(z)  given  by  equations  (106)-(103)  are  shown  in  Figures  21,  22 
and  23.  Comparison  of  these  results  shows  that  the  present  method  gives  better 

results . 
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Figure  21.  Nyquist  plot  of  the  frequency  responses  of  the  middle  loop 

of  the  continuous  system  and  the  corresponding  digital  control 
system  designed  by  different  methods  (example  2) . 
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Figure  23.  Phase  plot  of  the  frequency  responses  of  the  middle  loop  of  the 
continuous  system  and  the  corresponding  digital  control  system 
designed  by  different  methods  (example  2). 
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D^(z)  Controller: 

From  Figure  13,  the  transfer  function  of  the  continuous  model  for  the 
design  of  digital  controller  D^(z)  can  be  written  as 

*  s 

a  (s)  F,(s)C.(s)  * 

3V  7  R(s)  1+Gc4(s)F2(s)C3(s) 

.649S8-+30.5S7+559.9S6-K5534v104S5-K37  54x10SS4-K1883x106S3  . 

S9+62S8+187.2S7+.2359x105S6+.1711x106S5+.6447x106SA+.1373x10?S3 

e  n  / 

+.4981x10  S+. 3417xlO°S-627 .1 
+.1021x10?S2+.3466x106S-332.3 

(109) 


In  view  of  Figure  14  and  equations  (70)  and  (71) ,  we  can  write 


1v  s  an(s) 

xii+JYii  =  7Tf  *<jr) 

1+jYJ 


2  —IT 
s=j-  tan  y- 


and 


(110) 


*21+iY21  =  z(jY)“f 


0076w5-2 . 43w4+138 . 8w3+75 . 34w2+. 65w 


V'+M . 95w4+15 . 74w3-44 . 29w2- . 74w- . 2  7 


w=jy 


(111) 


where 


*(jY> 


D2(w) 


l+Gh0G1(w)D1(w)+Gh0G1G2(w)D2(w)D3(w) 


W=jy 


(112) 


Loading  equations  (109)-(112.)  into  the  computer  program  written  for 
the  algorithm  implementing  equations  (36)— (AO) ,  the  pulse-transfer  function 
of  the  digital  controller  D^Cz)  is  given  by 

D4(z)  «  1.0124  (113) 
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The  pulse-transfer  functions  of  D^(z)  obtained  by  the  existing  methods 
are  given  by 

Tustin  transform: 

D  (2)  «  .5858z2+.10652-.479 
4  z2-l. 082+0. 29 

Discrete  frequency  method: 

_  ,  v  1.994z2-3.716z+1.803  /11CX 

D  (z)  -  -r - -  ,  (115) 

*  z  -1.8812z+.9577 

The  frequency  responses  of  the  multi  loop  continuous  system  shown  in 
Figure  13  and  the  digital  control  system. shown  in  Figure  14  with  digital 
controller  D^(z)  given  by  equations  (113)-(115)  are  shown  in  Figures 
24-26.  Comparison  of  these  results  demonstrate  the  superiority  of  the 
present  method  over  the  existing  methods. 

SECTION  V 
CONCLUSIONS 

A  method  for  converting  an  existing  multiloop  continuous  control  system 
into  a  digital  control  system  with  similar  performances  is  developed.  An 
attempt  is  made  to  match  the  continuous  frequency  response  of  the  digital 
control  system  as  closely  as  possible  with  that  of  the  continuous  model.  The 
parameters  of  the  digital  controllers  are  obtained  by  minimizing  modified 
mean-square  of  the  errors  between  the  transfer  functions  of  the  continuous 
and  the  digital  systems  in  the  w-doraain.  Formulas  for  computing  the 
parameters  of  the  digital  controllers  are  obtained  and  have  been  programmed 
for  a  computer.  Since  the  design  procedure  is  computer-aided,  one  may  try  out 
several  digital  controllers  of  different  order  and  select  the  best  design  by 
considering  the  trade-off  between  the  cost  and  performance. 

The  numerical  examples  considered  show  that  the  performance  of  the 
digital  control  system  obtained  by  the  method  presented  is  superior  to  the 
ones  obtained  by  the  existing  methods. 
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Phase  plot  of  the  frequency  response  of  the  multiloop  continuous 
model  and  the  digital  control  system  designed  by  different 
methods  (example  2) . 


APPENDIX  A 


LISTING  OF  TRANSFER  FUNCTIONS 


In  this  appendix,  the  transfer  functions  of  the  aircraft  model 
(M  »  0.6,  h  =  30,000  ft)  used  in  Section  IV  are  provided. 


20 

s+20 


NUMERATOR 


I  NPOLY(I)  ZERO(I) 

1  (  -1.449  )S**  3  (  -.8003E-02)  +  J(  -.7455E-01) 

2  (  -147.6  )S**  2  (  -.8003E-02)  +  J(  .7455E-01) 

3  (  -2.370  )S**  1  (  -101.9  )  +  J(  0.  ) 

4  (  -.8295  ) 


DENOMINATOR 


I  DPOLY(I)  POLE (I) 

1  (  1.000  )S**  5  (  1.137  )  +  J(  0.  ) 

2  (  20.95  )S**  4  (  -.9345E-02)  +  J(  .7705E-01) 

3  (  16.71  )S**  3  (  -.9345E-02)  +  J(  -.7705E-01) 

4  (  -46.66  )S**  2  (  -2.071  )  +  J(  0.  ) 

5  (  -.7820  )S**  1  (  -20.00  )  +  J(  0.  ) 

6  (  -.2837  ) 


s+20 


NUMERATOR 


I  NPOLY(I)  ZERO (I) 

1  (  -147.6  )S**  3  (0.  )  +  J(  0.  ) 

2  (  -79.39  )S**  2  (  -.8816E-02)  +  J(  0.  ) 

3  (  -.6884  )S**  1  (  -.5292  )  +  J(  0.  ) 

4(0.  ) 


DENOMINATOR 


I  DPOLY(I)  POLE 


1 

(  1.000 

)S**  5 

(  .  1.137  )  +  J( 

0. 

) 

2 

(  20.95 

)S**  4 

(  — .-9345E-02)  +  J( 

.7705E-01) 

3 

(  16.71 

)S**  3 

(  - .9345E-02)  +  J( 

- .7705E-01) 

4 

(  -46.66 

)S**  2 

(  -2.071  )  +  J( 

0. 

) 

5 

6 

(  -.7820 
(  -.2837 

)S**  1 
) 

(  -20.00  )  +  J( 

0. 

) 
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